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Uber den funktionalen Zusammenhang zwischen 
Populationsdichte, Ausbreitungsvermögen und Fangmenge 
bei Bodenfallen 


Von Hubert Fechter 


Zoologische Staatssammlung Miinchen 


Abstract 


On the functional relation between population density, spreading activity and catch in 
pitfall traps. 

Starting with a partial differential equation which describes the two-dimensional, time- 
dependent density distribution of objects performing a random walk, the quantitative re- 
lation between population density, spreading activity, trap aperture and catch per time 
unit, is derived under the assumption of symmetrical isotropic and anisotropic conditions 
and the boundary condition of constant population density at the trap’s margin. 

The diffusion coefficient is introduced as a measure of the spreading activity and possi- 
bilities to determine it are considered. 

In order to keep the mutual influence and the reduction of the population density within 
certain limits when several traps are used, a formula is given which can be used to evaluate 
the reasonable distance of traps in a grid arrangement. 

The knowledge of the quantitative relation between the parameters is applicable to the 
determination of population density, of the spreading activity, the estimation of the 
capture chances of animals waiting for prey to approach and of the rate of immigration 
and transmigration in circular areas, as well as to some aspectes of island colonization. 

To simplify the exploitation of the functional relation between the variables a table of 
associated values is provided. 


Einleitung 


Populationsdichte und Ausbreitungsvermögen sind Größen, deren Kenntnis bei 
vielen ökologischen Untersuchungen, insbesondere bei der Analyse populationsdy- 
namischer Vorgänge,.von entscheidender Bedeutung ist. Zu ihrer Feststellung sind 
viele Methoden erdacht worden, die alle auf Fangtechniken basieren, deren techni- 
sche Einzelheiten wesentlich von der Lebensweise der zu untersuchenden Tiere be- 
stimmt werden. Zum Fang von auf dem Boden umherlaufenden, flugunfähigen Tie- 
ren, wie zahlreichen Käfern, Tausendfüßlern, Spinnen und Asseln werden fast aus- 
schließlich Bodenfallen — nach BARBER (1931), der sie erstmals systematisch ein- 
setzte, auch Barber-Fallen genannt — verwendet. Es sind dies fallgrubenartig in 
den Boden eingesenkte, in der Regel zylindrische Gefäße mit oder ohne Köder be- 
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ziehungsweise Konservierungsflüssigkeit. Tiere, die beim Umherlaufen über den 
Rand der Bodenfalle geraten stürzen hinein und können bei richtig konstruierter 
Falle nicht mehr entkommen. 

Die Zahl der gefangenen Tiere hängt von mehreren Faktoren ab: 1. den Ausma- 
ßen der Falle, 2. der Populationsdichte, 3. dem Ausbreitungsvermögen und der Ak- 
tivität der Tiere sowie 4. dem Wirkungsgrad der Falle. Ihr Einfluß auf das Fang- 
ergebnis und die daraus möglichen Rückschlüsse wurden bereits häufig untersucht 
und diskutiert (TRETZEL 1955, HEYDEMANN 1953, 1958, 1962; MITCHELL 1963, 
GREENSLADE 1964, Hinps & Ricxarp 1973, Lurr 1973, 1975 u. a.), jedoch erfolgte 
bisher keine exakte Analyse der theoretischen Beziehungen zwischen den oben ge- 
nannten Parametern. Vielen erscheinen diese Zusammenhänge als logisch zwingend 
gegeben, andere dagegen (HEYDEMANN 1958, BRIGGS 1961) bestreiten sie. 

Ziel dieser Arbeit ist es die Abhängigkeit der Zahl der Fallenfänge pro Zeitein- 
heit von der Populationsdichte, dem Ausbreitungsvermögen und der Offnungsweite 
köderloser Fallen in ihrem quantitativen Zusammenhang zu untersuchen und ihre 
möglichen praktischen Nutzanwendungen zu diskutieren. 


Theoretische Zusammenhänge 


Die zeitliche Abhängigkeit der Konzentrationsverteilung von Objekten, die eine ` 
Zufallsbewegung (random walk) ausführen, läßt sich, wie CHANDRASEKHAR (1943) 
gezeigt hat, durch eine partielle Differentialgleichung beschreiben, die, für die Ebene 
und radialsymmetrische Bedingungen im Zylinderkoordinaten ausgedrückt, folgen- 
dermaßen lautet 
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wobei C die Konzentration, d. h. im vorliegenden Falle die Flächendichte der Indi- 
viduen, r den radialen Abstand vom Koordinatenursprung — hier zweckmäßiger- 
weise die Entfernung von der Achse des Fallenhohlraumes im Niveau der Boden- 
oberfläche —, t die Zeit und D ein Maß für das Ausbreitungsvermögen bedeuten. 

Die Lösung von (1) wird wesentlich von den herrschenden Randbedingungen be- 
stimmt. Die gewählten Randbedingungen stellen einen annehmbaren Kompromiß 
zwischen den biologischen Gegebenheiten und der zur mathematischen Formulie- 
rung und Behandlung notwendigen Abstraktion und Idealisierung dar. Die Popula- 
tionsdichte wird als hinreichend groß vorausgesetzt, derart, daß durch den Entzug 
der gefangenen Tiere der Dichtewert während der Fangperiode in der unmittelba- 
ren Umgebung des Fallenrandes nicht wesentlich verändert wird, die Populations- 
dichte am Fallenrand also im Mittel als konstant betrachtet werden kann. Die Popu- 
lation soll ferner eine endliche Ausdehnung besitzen, d. h. ihre Dichte in unendlicher 
Entfernung von der Falle gleich null sein. Der Fallenrand bei r =a wird im Sinne 
des random walk als absorbierende Barriere betrachtet. Zusammengefaßt lauten die 
Randbedingungen also: 


1. Beit > O undr = aist C = C, = const. 
2. Beit >Oundr = œ ist C = O. 


Die Lösung der partiellen Differentialgleichung (1) erfolgt zweckmäßigerweise 
analog zu der formal gleichen Differentialgleichung der Wärmeleitung (Taurz 
1971) mit Hilfe der LapLace-Transformation 
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und ihrer Umkehrung. Durch die Transformation von (1) ergibt sich 
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eine gewöhnliche Differentialgleichung. Die transformierte Randbedingung C = C, 
lautet u = C,/p. 
Ein Lösungsansatz, der die Bedingungen erfüllt, ist 
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wobei H,(') eine modifizierte BesseL-Funktion, die sogenannte HaAnker’sche 
Funktion O. Ordnung darstellt, die als einzige mit unendlich werdendem Argument 
gegen null geht. 

Die Integrationskonstante A ergibt sich aus der transformierten Randbedingung 
u= C,/p beir = a, dem Radius des Fanggefäßes, zu 
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und die vollständige Lösung von (2) lautet: 
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Da die Lösungsfunktion (4) bei p = O einen Verzweigungspunkt besitzt, muß 
die Rücktransformation mit Hilfe des Residuensatzes über einen geschlossenen Weg 
erfolgen, der die Singularitäten umgeht und ausschließt. Der Integrationsweg ver- 
läuft dabei entlang einer Parallelen im Abstand a zur imaginären Achse so, daß 
die singuläre Stelle des Integranden links von ihr liegt, geht über den oberen Vier- 
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telkreis der linken Halbebene zum oberen Ufer eines entlang der negativen p-Achse 
gelegten Verzweigungsschnittes, läuft an diesem entlang zum Nullpunkt, der in 
einem kleinen Kreis umgangen wird und führt über das untere Schnittufer und den 
unteren Viertelkreis zur Geraden mit dem Abszissenwert a zurück. Die Integrale 
über die Viertelkreise verschwinden, wenn deren Radius gegen Unendlich geht. Da 


die Summe der Integrale über die Berandung gleich null sein muß, a+i co 
kann das Bee 


durch die Summe der Integrale über die Wege entlang des positiven und negativen 
Schnittufers und des Kreises um den Ursprung ersetzt werden. 
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Das Umlaufintegral $ um den Nullpunkt ergibt den Wert 1, wenn der Radius des 


Umlaufkreises gegen O geht, da sowohl der Quotient aus den Zylinderfunktionen 
als auch die Exponentialfunktion diesen Grenzwert annehmen und auch 


Setzt man p/D1 = q, und da der Verzweigungsschnitt entlang der negativen 
p-Achse verläuft p = -Dq?, und berücksichtigt, daß am oberen Schnittufer das po- 
sitive, am unteren das negative Vorzeichen gilt, sowie dp/p = 2dq/q ist, so ergibt 
(5) auf (4) angewandt 


oo (0) 
een Oe) a, ee 1 | HC) han) 
Se em ae oe er | gC") (Fea) 
e-Da’t 2dq l (6) 
q 


Zwischen den Hanxetschen Funktionen und den Besset- [ J (qr) | - 
beziehungsweise Neumann- [ N, (qr) | -Funktionen gelten folgende Beziehungen: 


Ea) (Iqrl) =) (qr) T iNo (qr) =) (-qr) + iNo (-qr) == [ Jo (qr) -iNo (qr) ] 


Führt man diese in (6) ein, so erhält man schließlich 


Tan (qa) N (gr) — Jo (ar) No (qa) — Jo (ar) No (qa) ek 
C= eee EZ 
Jo (aa) + N3 (qa) q | r 


die Funktion, aus der a die Fangmenge berechnen läßt. 
Die Anzahl der sich zufällig bewegenden Tiere, die in der Zeit t bei r = a den 
Fallenrand von der Lange s = 2 ra überschreiten beträgt 
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Die Differentiation von (7) nach r unter Beriicksichtigung der Beziehungen 
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a lp (91) = -J, (ar) bzw. en N, (qr) = -Ni(qr) und 
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sowie die Multiplikation entsprechend (8) ergibt. 
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Das in (9) auftretende Integral kann numerisch ausgewertet werden. Für große 
Funktionswerte, wie sie im biologischen Bereich vorkommen, hat JAEGER (1942) 
eine Formel zur näherungsweisen Berechnung des Integrals angegeben, die, in (9) 
eingesetzt, zu folgendem Endergebnis führt: 


1 0,57722 
re nn nn EP 
[In (4Dr/a2)- 1,15444] [In(4Dt/a2)-1,15444] 


s 1,31175 A 1 k 
[In (4Dr/a?) - 1,15444 J? [In (4Dt/a2) - 1,15444 ]* (10) 


Der Ausdruck in den geschweiften Klammern ist für Werte von Dt/a? zwi- 
schen 10? und 105 in Tafel I tabelliert. 

Wie aus (10) hervorgeht, besteht zwischen der Fangmenge M, gemessen in Indi- 

viduen, dem Offnungsradius a der Falle in Metern, der Populationsdichte C, [In- 
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dividuen/m?] und dem Ausbreitungsvermögen D [m?/Tag] eine strenge Abhängig- 
keit. Von all diesen Größen bedarf das Ausbreitungsvermögen einer näheren Er- 
läuterung. 


Tafel 1: Tabelle der Werte des Ausdrucks in t der Gleichung (10) für Dt/a? zwischen 
10? und 105.) 
00 10 20 30 40 50 60 70 80 90 


0,1723 0,1698 0,1675 0,1655 0,1637 0,1620 0,1604 0,1590 0,1577 0,1564 
0,1553 0,1542 0,1532 0,1522 0,1513 0,1504 0,1496 0,1488 0,1480 0,1473 
0,1466 0,1460 0,1453 0,1447 0,1441 0,1436 0,1430 0,1425 0,1420 0,1415 
0,1410 0,1405 0,1401 0,1397 0,1392 0,1388 0,1384 0,1380 0,1377 0,1373 
0,1369 0,1366 0,1362 0,1359 0,1356 0,1352 0,1349 0,1346 0,1343 0,1340 
0,1337 0,1335 0,1332 0,1329 0,1326 0,1324 0,1321 0,1319 0,1316 0,1314 
0,1312 0,1309 0,1307 0,1305 0,1302 0,1300 0,1298 0,1296 0,1294 0,1292 
0,1290 0,1288 0,1286 0,1284 0,1282 0,1280 0,1279 0,1277 0,1275 0,1273 
0,1271 0,1270 0,1268 0,1266 0,1265 0,1263 0,1261 0,1260 0,1258 0,1257 


O oo NGAUN È L N e 


000 100 200 300 400 500 600 700 800 900 


0,1255 0,1241 0,1228 0,1217 0,1206 0,1197 0,1188 0,1180 0,1172 0,1165 
0,1158 0,1152 0,1146 0,1140 0,1135 0,1130 0,1125 0,1120 0,1116 0,1112 
0,1108 0,1104 0,1100 0,1096 0,1093 0,1089 0,1086 0,1083 0,1080 0,1077 
0,1074 0,1072 0,1069 0,1066 0,1064 0,1061 0,1059 0,1056 0,1054 0,1052 
0,1050 0,1048 0,1046 0,1044 0,1042 0,1040 0,1038 0,1036 0,1034 0,1032 
0,1031 0,1029 0,1027 0,1025 0,1024 0,1022 0,1021 0,1019 0,1018 0,1016 
0,1015 0,1013 0,1012 0,1011 0,1009 0,1008 0,1007 0,1005 0,1004 0,1003 
0,1001 0,1000 0,0999 0,0998 0,0997 0,0996 0,0994 0,0993 0,0992 0,0991 
0,0990 0,0989 0,0988 0,0987 0,0986 0,0985 0,0984 0,0983 0,0982 0,0981 


oc. ON BD wT A L Ne 


0000 1000 2000 3000 4000 5000 6000 7000 8000 9000 


0,09800 0,09711 0,09632 0,09559 0,09493 0,09433 0,09377 0,09325 0,09276 0,09231 
0,09188 0,09148 0,09110 0,09074 0,09039 0,09007 0,08976 0,08946 0,08917 0,08890 
0,08864 0,08838 0,08814 0,08791 0,08768 0,08746 0,08725 0,08704 0,08685 0,08665 
0,08647 0,08629 0,08611 0,08594 0,08577 0,08561 0,08545 0,08530 0,08514 0,08500 
0,08485 0,08471 0,08458 0,08444 0,08431 0,08418 0,08406 0,08393 0,08381 0,08370 
0,08358 0,08347 0,08335 0,08325 0,08314 0,08303 0,08293 0,08283 0,08273 0,08263 
0,08253 0,08244 0,08234 0,08225 0,08216 0,08207 0,08198 0,08190 0,08181 0,08173 
0,08164 0,08156 0,08148 0,08140 0,08133 0,08125 0,08117 0,08110 0,08102 0,08095 
0,08088 0,08081 0,08074 0,08067 0,08060 0,08053 0,08046 0,08040 0,08033 0,08027 
0,08020 
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*) Beispiel: Für Dt/a? = 3400 ist n = 0,1093 


Das Ausbreitungsvermögen 


Für ein Tier, das sich kaum fortbewegt, ist die Wahrscheinlichkeit auf die Falle zu 
stoßen ziemlich gering, für ein sehr vagiles Tier dagegen ist die Gefahr in die Falle 
zu laufen erheblich größer. Dabei kommt es nicht darauf an, welche Entfernung das 
Tier zurücklegt, sondern wie groß das Gebiet ist, in dem es sich während einer be- 
stimmten Zeit theoretisch bewegen kann. Gesucht ist daher ein Maß für die poten- 
tielle, flächenbeherrschende Fortbewegungsaktivität, mit einem Wort für das Aus- 
breitungsvermögen des Tieres. Aus der Theorie des Random Walk und der Brown- 
schen Bewegung ergibt sich, ganz unabhängig von der Natur der sich bewegenden 
Objekte, als Maß hierfür der Diffusionskoeffizient, der häufig auch Diffusionskon- 
stante oder Diffusivität, gelegentlich auch Dispersität genannt wird. 

Bei rotationssymmetrischen Bedingungen und für die Ebene betrachtet ist der Dif- 
fusionskoeffizent D eines Individuums dem mittleren Verschiebungsquadrat L? 
während eines definierten Zeitabschnittes (t) proportional, wobei sich aus statisti- 
schen Überlegungen ein Proportionalitätsfaktor von !/4 ergibt (PrELou 1969). 


Do Zu WE (11) 


L? selbst ist das Mittel aus der Summe der Quadrate der Längen an der Ver- 
bindungsgeraden zwischen den jeweiligen Anfangs- und Endpunkten der Bewegung 
innerhalb eines gewählten, konstanten Zeitintervalls (r<t); bein gemessenen Strek- 
ken ist demnach 
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Voraussetzung dabei ist, daß das gewählte Zeitintervall (t) im Verhältnis zum 
gesamten Beobachtungszeitabschnitt (t = n t)sehr klein, d. h. n sehr groß ist und 
bei den von dem Objekt ausgeführten Zick-Zack-Bewegungen alle Richtungen 
gleich häufig eingeschlagen werden. Unter diesen Bedingungen ist D isotrop — in al- 
len Richtungen der Ebene gleich groß. 

Die Belaufareale der Tiere einer Population sind jedoch nur selten kreisförmig, 
meist sind die Aktionsräume (SCHWERDTFEGER 1968), den Substratbedingungen ent- 
sprechend, unregelmäßig begrenzt, lassen sich aber häufig durch eine elliptische Kon- 
figuration annähern. In diesem Falle ist D anisotrop, d. h. in Richtung der beiden 
Ellipsenachsen gesehen unterschiedlich groß. Die beiden Diffusionskoeffizienten D, 
und D, lassen sich dann analog zu dem oben geschilderten Verfahren aus den 
Quadraten der Projektionen der an auf die X- beziehungsweise Y-Achse ermitteln, 
wobei man es mit zwei eindimensionalen Random Walks auf den Achsen mit unter- 
schiedlichen mittleren Verschiebungslängen zu tun hat, für die dann in der Bezie- 
hung gemäß (11) zwischen D, und Li respektive D, und L ein Proportiona- 
litatsfaktor von 1/2 gilt. 

Die Differentialgleichung (1) fiir den anisotropen Fall lautet dann 
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Durch Transformation auf die neuen Variablen & = x ]/D/D, und 
N | D/D, läßt sich (13) jedoch leicht in den isotropen Fall 
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überführen, so daß sich am Endergebnis (10) prinzipiell nichts ändert, außer daß der 
Diffusionskoeffizient durch D = }]/ D, Dy ersetzt werden muß. 


Im Ausbreitungsprozeß läßt sich die Wahrscheinlichkeit, mit der sich das Objekt 
an einem ganz bestimmten Ort aufhält, mit Hilfe einer Normalverteilung, in der 
das Produkt aus Diffusionskoeffizient und Zeit (2Dt) als Varianz auftritt, darstel- 
len (PreLou 1969). Dies liefert auch einen Ansatzpunkt zur experimentellen Bestim- 
mung des Diffusionskoeffizienten in einer Population. Nach den von SKELLAM 
(1951) angestellten Berechnungen über die Ausbreitungsgeschwindigkeit von Popu- 
lationen besteht zwischen dem Radius des Grenzkreises R,, bis zu dem sich nach 
Ablauf der Zeit t die an einem bestimmten Ausgangspunkt freigelassene Anzahl von 
N Tieren ausgebreitet hat, und dem Diffusionskoeffizienten die Beziehung R? = 
4DtInN; dabei darf das am weitesten vorgedrungene Individuum nicht mitgerech- 
net werden. Hieraus folgt: 


Re 


a 4t In N 


(14) 


Diese Methode wäre auch auf den anisotropen Fall anwendbar, da bei einer ellip- 
tischen Ausbreitung die Haupt- und Nebenkreis-Radien R,. und R, , der Ellipse 
zur Ermittlung von D, beziehungsweise D, nach (14) herangezogen werden 
könnten und 
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den in (10) einsetzbaren Diffusionskoeffizienten ergäbe. 

Andere Verfahren zur Bestimmung achsenbezogener Diffusionskoeffizienten bei 
elliptischer Ausbreitung finden sich bei CLARK (1962) und Paris (1965). Eine wei- 
tere Möglichkeit bietet das Capture-Recapture-Verfahren (BAILEY 1951) mit sei- 
nen Abwandlungen (Jolly-Methode und Bailey’s triple catch; vgl. MÜHLENBERG 
1976) zur Feststellung der Populationsgröße in Arealen, die gut abgrenzbar sind 
und in denen sich durch die Anordnung der Fallen eine akzeptable Zuordnung von 
Individuenzahlen und Größe der Siedlungsfläche herstellen und dadurch die Popu- 
lationsdichte feststellen läßt. Der Diffusionskoeffizient könnte dann aus (10) durch 
iterative Verfahren, mit Hilfe von Tabellen (siehe Tafel I) oder Nomogr aune] 
gewonnen werden. 


Anordnung der Fallen 


Zur Steigerung der Ausbeute und um ein Gebiet möglichst gleichmäßig zu befan- 
gen, ist es üblich, gleichzeitig mehrere Bodenfallen, meist in gitterförmiger Anord- 
nung, zu verwenden. Nur um Tiere zu erbeuten, sich einen Überblick über die Fau- 
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na zu verschaffen oder für das Capture-Recapture-Verfahren ist dies durchaus 
zweckmäßig. Soll jedoch mit Hilfe der Fänge einer der Parameter in (10) ermittelt 
werden, so ist die gegenseitige Beeinflussung der Fallen zu berücksichtigen. Wie 
Lurr (1975) gezeigt hat, verdecken sich die Fallen gegenseitig in bestimmten Sek- 
toren, wodurch ihre Effektivität herabgesetzt wird. Von den in (10) auftretenden 
Parametern können die Populationsdichte und der Diffusionskoeffizient durch An- 
zahl und Anordnung der Fallen verändert werden. Die wesentliche, eingangs aufge- 
stellte Randbedingung war die Konstanz der Populationsdichte während der Dauer 
der Fangperiode. Das bedeutet, die Zahl der Fallen muß so gering und die Fang- 
periode so kurz gehalten werden, daß durch den Entzug der gefangenen Tiere die 
Populationsdichte am Fallenrand im Mittel nicht wesentlich verringert wird und die 
benachbarten Fallen sich gegenseitig keine Konkurrenz machen. Da sich die beiden 
Störgrößen — Verringerung der Populationsdichte durch Entzug der gefangenen 
Tiere der untersuchten Art und die Konkurrenz der Fallen — nie vollständig aus- 
schalten lassen, gilt es, sie so klein zu halten, wie die Versuchsbedingungen es erfor- 
dern. 

Die Wahrscheinlichkeit, daß ein Individuum, das am Fallenrand startet, in der 
Zeit t über einen Grenzkreis vom Radius R, hinaus gelangt, beträgt (vgl. PıELou 
1969) 


Ri 
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Pr = exp G 


Sind die Fallen an den Knotenpunkten eines R,X R, Gitternetzes angeordnet, 
so daß auf den Grenzkreisen im 90°-Abstand jeweils eine Falle mit dem Off- 
nungsradius a liegt, so ist die Wahrscheinlichkeit dafür, daß ein den Grenzkreis 
überschreitendes Tier dabei auf eine der 4 Fallen stößt 
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TR: 


Pe = (16) 


Für das gleichzeitige Eintreten beider Ereignisse, d. h. Grenzüberschreitung und 
Fallenkontakt, gilt: 


a R; 
ap e=) 17 
TRe p 4Dt ae 


Pr = h S 


Mit pp haben wir ein Maß für die Unsicherheit, die aus der Konkurrenz der 
Fallen untereinander erwächst. 

Der Einfluß des Entzugs der gefangenen Tiere läßt sich folgendermaßen abschät- 
zen: Ist M das Mittel der während der Fangperiode pro Falle gefangenen Indivi- 
duen einer Art und F die Fläche mit den Seitenlängen g und h auf der N Fallen auf- 
gestellt sind, so ergibt sich der Dichteverlust der Population zu 
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NG = Gye oz (18) 


Auf der Fläche F = gh lassen sich in einer R, X R, Gitteranordnung 
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5 
= (— - = 19 
N= R, RO R, R.) (19) 
Fallen unterbringen; (19) in (18) führt auf 
h 
AC = ne) ug = Ry) (— = R.) (20) 
gh Re t 


Es sei nun die Wahrscheinlichkeit dafür, daß eine bestimmte Anzahl von Indi- 
viduen pro Flächeneinheit entzogen wird gleich py. , d.h. AC = Pac Co; aus (20) 
folgt dann: 
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Setzt man die Wahrscheinlichkeit P = pae Pk, mit der beide Störtaktoren gleich- 
zeitig vorkommen dürfen, von vornherein fest, so läßt sich der unter den gegebenen 
Bedingungen einzuhaltende Fallenabstand R, aus 


M g h 4a R? 
P = 12 en (Se X = 22 
| Cash x ) Re | RR: N 7D: ? = 
iterativ berechnen, beziehungsweise es läßt sich nachträglich ermitteln, welcher 
Unsicherheitsfaktor bei bestimmten Gegebenheiten vorgelegen hat. 


Diskussion und Anwendungsmöglichkeiten 


Die Kenntnis der genauen quantitativen Beziehungen erlaubt es, im Gegensatz zu 
den bisherigen Verfahren, erstmals Berechnungen der abhängigen Variablen vorzu- 
nehmen und in Kombination mit anderen Methoden die experimentelle Bestimmung 
populationsdynamisch wichtiger Parameter durchzuführen. 

Geht man von den Gleichungen (10) beziehungsweise (22) aus, so haben wir es 
mit fünf, respektive acht Variablen zu tun, von denen im ersten Fall zwei, nämlich 
die Zeit und der Fallenöffnungsradius a, im zweiten Falle vier, das sind die Areal- 
ausmaße g, h, die Wahrscheinlichkeit P und ebenfalls die Zeit t, in gewissen Gren- 
zen frei gewählt und dann konstant gehalten werden können. Übrig bleiben drei 
Variable: Fangmenge M, Populationsdichte C, und Diffusionskoeffizient D — be- 
zieht man das Fallenancrdnungsproblem mit dem Fallenabstand R, ein, so kommt 
noch eine vierte hinzu —, von denen jeweils eine als abhängige Veränderliche be- 
trachtet und mit Hilfe der restlichen berechnet werden kann, wenn diese bekannt 
sind oder in einer bestimmten Größenordnung angenommen oder vorausgesetzt 
werden. 

Am relativ unkritischsten dürfte der Fallenabstand R, sein, für dessen Ermitt- 
lung grobe Schätzwerte der anderen Parameter, die sich für diesen Zweck in vielen 
Fällen aus der Literatur ableiten lassen, als ausreichend angesehen werden könnten. 

Da die Fangmenge als Versuchsergebnis stets gegeben ist, haben wir es vor allem 
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mit den Unbekannten D und C, zu tun. Ist D bekannt, so läßt sich aus dem Fang- 
ergebnis direkt die Populationsdichte errechnen und umgekehrt bei bekanntem C, 
der Diffusionskoeffizient. Methoden zur Bestimmung beider Größen wurden oben 
bereits erörtert. 

Um zu den tatsächlich verrechenbaren Fangzahlen zu kommen, ist die Fangmen- 
ge M, wie sie sich aus (10) ergibt noch mit einem Korrekturfaktor, der dem rezipro- 
ken Wert des Nutzeffektes der Falle entspricht, zu multiplizieren. In jüngster Zeit 
hat Lurr (1975) den Wirkungsgrad von Bodenfallen bezüglich ihrer Fang- und 
Rückhalteeffizienz ın Abhängigkeit von den Fallenausmaßen, der Materialbeschaf- 
fenheit und der Tiergröße untersucht. Seine Versuche zeigen, daß bei Käfern nur et- 
wa 51°/o bis 879/o der Tiere, die mit dem Fallenrand in Kontakt kommen auch ge- 
fangen werden; hinsichtlich der Fangeffizienz scheint zwischen der Fallenöffnungs- 
weite und der Körpergröße der Tiere ein umgekehrt proportionales Verhältnis zu 
bestehen. Es sollte deshalb in Vorversuchen ein zweckmäßiger Kompromiß zwi- 
schen dem Durchmesser der Falle und ihrer Fangeffizienz gefunden werden. Die 
Entkommensrate ist bei Glaswänden so gering, daß sie praktisch vernachlässigt wer- 
den kann. 

Die Anwendungsmöglichkeiten, die die Kenntnis der quantitativen Zusammen- 
hänge bietet, beschränken sich nicht nur auf die Fangtechnik mit Bodenfallen. Es 
lassen sich damit unter anderem die Fangchancen von Lauerern, beispielsweise Fall- 
türspinnen, und die für ihre ausreichende Ernährung notwendige Populationsdichte 
der Beutetiere abschätzen oder die Einwanderungsrate in ein bestimmtes kreisför- 
ıniges Areal, z. B. auch Inseln, beziehungsweise die Anzahl der Durchzügler ange- 
ben. 

Das Hauptanwendungsgebiet sehe ich jedoch in der Bestimmung des Diffusions- 
koeffizienten in Kombination mit dem Capture-Recapture-Verfahren zur Feststel- 
lung der Populationsdichte. Es gibt zwar eine Fülle von Untersuchungen über Ak- 
tivitätsphasen und Laufgeschwindigkeiten von Bodentieren, bisher aber keinerlei 
experimentelle Bestimmung der sehr viel aussagekräftigeren Diffusionskoeffizien- 
ten, deren Kenntnis bei der Verifizierung vieler populationsdynamischer Modelle 
wie Migrationstheorien (KERNER 1959), genetischen Durchmischungsvorgängen und 
in Zukunft sicherlich auch bei der Bewältigung der im Natur- und Umweltschutz 
auftretenden Probleme und den in diesem Zusammenhang notwendigen modellhaf- 
ten Untersuchungen unentbehrlich ist. 

Der Diffusionskoeffizient einer Art ist eine zeitlich wie räumlich nicht konstante 
Größe, die ebenso wie die Aktivität beträchtlichen Schwankungen unterliegen kann 
und für die sich nur unter ganz bestimmten Bedingungen geltende Mittelwerte wer- 
den finden lassen. Diese Schwankungen von D sind letztlich auch die Ursache für 
den bei Fallenfängen häufig vorgetäuschten raschen Dichtewechsel in Populationen, 
der HEYDEMANN (1953, 1962) zur Prägung des Begriffes der Aktivitätsdichte ver- 
anlaßte. Faktoren, von denen ein Einfluß zu erwarten ist, sind Alter, Geschlecht, 
Entwicklungszustand, Bodenbeschaffenheit und Vegetationsbestand, Temperatur, 
Luftfeuchtigkeit, kurz Witterungseinflüsse, saisonale sowie circadiane Aktivitats- 
rhythmen und möglicherweise eine Dichteabhängigkeit. Einige der Einflüsse könn- 
ten sich, über die Fangperiode hin gesehen, gegenseitig aufheben oder herausmitteln 
und es ließen sich so wochen- oder tagescharakteristische D-Werte für definierte 
Biotope gewinnen. Die vorliegende, quantitative Analyse könnte wesentlich zur 
Aufklärung des Einflußgefüges beitragen. 
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Zusammenfassung 


Ausgehend von einer partiellen Differentialgleichung, die die zweidimensionale 
Dichteverteilung von Zufallsbewegungen (random walk) ausführenden Objekten in 
Abhängigkeit von der Zeit beschreibt, wird für Bodenfallen, unter der Annahme 
rotationssymmetrischer, aber auch radialsymmetrisch-anisotroper Bedingungen so- 
wie der Randbedingung konstant bleibender Dichte am Fallenrand, der quantitati- 
ve Zusammenhang zwischen Populationsdichte, Ausbreitungsvermögen, Offnungs- 
weite der Falle und Fangmenge pro Zeiteinheit abgeleitet. 

Der Diffusionskoeffizient als Maß des Ausbreitungsvermögens wird erläutert und 
experimentielle Bestimmungsmöglichkeiten dafür aufgezeigt. 

Um bei gleichzeitiger Verwendung mehrerer Fallen die gegenseitige Beeinflus- 
sung und Verminderung der Populationsdichte in gewünschten Grenzen halten zu 
können, wird ein Verfahren zur Berechnung des zweckmäßigen Fallenabstandes 
angegeben. 

Die Kenntnis des quantitativen Zusammenhanges der Parameter bietet die Mög- 
lichkeit zur Bestimmung von Populationsdichten, vor allem aber des Ausbreitungs- 
vermögens, gestattet das Abschätzen der Fangchancen von Lauerern sowie der Ein- 
und Durchwanderungsraten in kreisförmigen Arealen und der Besiedlungsvorgänge 
von Inseln. 

Eine Tabelle mit korrespondierenden Zahlenwerten soll die Auswertung der 
funktionalen Beziehungen zwischen den Veränderlichen vereinfachen. 
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